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•  A  novel  normalized  OCV  modeling  approach  is  developed. 

•  The  normalized  OCV  modeling  is  stable  over  temperature  changes  and  battery  aging. 

•  The  importance  of  battery  capacity  in  OCV  modeling  is  demonstrated. 

•  The  hysteresis  is  marginalized  by  having  symmetric  charging  and  discharging  data. 

•  Proposed  OCV  modeling  is  validated  for  different  chemistries,  temperatures  and  aging. 
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Several  aspects  of  the  open  circuit  voltage  (OCV)  characterization  of  Li-ion  batteries  as  it  applies  to 
battery  fuel  gauging  (BFG)  in  portable  applications  are  considered  in  this  paper.  Accurate  knowledge  of 
the  nonlinear  relationship  between  the  OCV  and  the  state  of  charge  (SOC)  is  required  for  adaptive  SOC 
tracking  during  battery  usage.  BFG  in  portable  applications  requires  this  OCV-SOC  characterization  to 
meet  additional  constraints:  (i)  The  OCV-SOC  characterization  has  to  be  defined  with  a  minimum 
number  of  parameters;  (ii)  It  should  be  easily  computable  and  invertible  with  few  operations;  and  (iii) 
Computation  of  the  model,  its  derivative  and  its  inverse  should  be  possible  in  a  numerically  stable  way. 
With  the  help  of  OCV-SOC  characterization  data  collected  from  34  battery  cells  each  at  16  different 
temperatures  ranging  from  -25  °C  to  50  °C,  we  present  the  following  results  in  this  paper:  (a)  A  robust 
normalized  OCV  modeling  approach  that  dramatically  reduces  the  number  of  OCV-SOC  parameters  and 
as  a  result  simplifies  and  generalizes  the  BFG  across  temperatures  and  aging,  (b)  Several  novel  functions 
for  OCV  modeling,  (c)  Efficient  methods  to  simplify  the  computations  of  OCV  functions,  (d)  Novel 
methods  for  OCV  parameter  estimation,  and  (e)  A  detailed  performance  analysis. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Rechargeable  Lithium-ion  batteries  are  widely  used  in  electric 
vehicles,  personal/wearable  electronic  devices,  household  appli¬ 
ances,  aerospace  equipment,  grid  storage,  and  so  on.  Estimating  the 
state  of  charge  (SOC)  and  state  of  health  (SOH),  known  as  battery 
fuel  gauging  [1—3],  is  an  important  function  in  equipment  powered 
by  rechargeable  batteries.  The  OCV  is  a  measure  of  the  electro¬ 
motive  force  (EMF)  of  the  battery,  which  is  known  to  have  a 
monotonic  relationship  with  the  SOC  of  the  battery,  hence,  esti¬ 
mating  SOC  must  have  been  a  straightforward,  voltage-look-up 
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process.  However,  due  to  the  variable  internal  resistance,  capac¬ 
ity,  hysteresis  and  relaxation  effects  of  the  battery,  the  BFG  involves 
many  components  as  outlined  in  Fig.  1;  please  see  Ref.  [2  for  a 
detailed  description  of  the  function  of  each  block.  The  focus  of  the 
paper  is  on  estimating  the  parameters  of  the  OCV-SOC 
characterization. 

Existing  OCV  modeling  approaches  can  be  broadly  classified  into 
chemistry-based  and  Current-Voltage  based  approaches.  In 
chemistry-based  approaches  [4  ,  the  OCV  of  each  electrode  (anode 
and  cathode  w.r.t.  some  reference)  is  expressed  as  a  function  of  the 
utilization  of  the  electrode  (the  lithium  concentration  in  the  elec¬ 
trode  normalized  by  the  maximum  possible  concentration)  or  the 
SOC  of  each  electrode.  It  is  generally  assumed  that  this  anode  and 
cathode  SOC  varies  linearly  with  the  cell  SOC  [5].  Subsequently,  the 
difference  between  the  OCV  of  the  anode  and  cathode  gives  the 
OCV  of  the  complete  cell.  High  current  rates  (i.e.,  near  the  rated 
maximum)  have  been  shown  to  affect  the  macroscopic  processes  in 
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Fig.  1.  Elements  of  a  BFG.  The  focus  of  this  paper,  OCV-SOC  characterization  is  shaded  in  gray.  Details  of  System  ID  block  can  be  found  in  Refs.  [1,2]  and  the  details  of  SOC  tracking 
block  can  be  found  in  Ref.  3]. 


a  way  that  the  OCV  hysteresis  vanishes  for  Li-ion  cells,  which 
regularly  show  OCV  hysteresis  after  low  current  application  [6]. 
Roscher  et  al.  [7  conducted  OCV  (full  and  partial  charge-discharge 
cycle)  tests  on  Lithium  ion  phosphate  (LiFePO^  batteries  to  char¬ 
acterize  the  hysteresis  and  recovery  effects.  The  final  OCV  model  is 
constructed  by  concatenating  the  actual  SOC,  the  recovery  factor 
and  the  hysteresis  factor.  Our  focus  in  this  paper  is  on  Cur¬ 
rent-Voltage  based  OCV-SOC  characterization. 

The  Current- Voltage  based  OCV-SOC  characterization  can  be 
summed  by  two  simple  steps: 

(i)  Collect  pairs  of  {OCV,  SOC}  values,  spanning  the  entire  range 
of  SOC=[0,1] 

(ii)  Use  the  above  data  to  estimate  the  parameters  of  the  func¬ 
tion  OCV  =j[ SOC)  for  a  hypothesized  function/. 

There  are  challenges  in  both  of  the  above  steps.  The  objective  of 
this  paper  is  to  detail  them  and  discuss  ways  to  address  them. 

Fig.  2  shows  the  equivalent  circuits  of  a  battery  when  it  is 
experiencing  dynamic  current  versus  when  it  is  rested  [1  .  When 
the  battery  is  experiencing  dynamic  current,  the  voltage  difference 
between  OCV,  V0(s[/c]),  where,  s[k\  is  the  SOC  at  the  discrete  time  k , 
and  the  measured  terminal  voltage  v[/<]  becomes  higher.  When  the 
current  i[k]  becomes  zero,  the  battery  starts  to  “relax”,  i.e.,  the 
capacitors  Ci,...,Qv  start  to  discharge  through  Rp  ...,  R !jy,  respec¬ 
tively.  The  battery  is  said  to  be  in  a  “rest  state”  when  i[k]  =  0  for 
longer  than  the  multiple  (typically  >5)  of  the  largest  time  constant 
RiQ  in  the  equivalent  circuit  model,  i.e.,  until  all  the  capacitors 
Ci,...,Qv  are  fully  discharged.  However,  the  OCV  is  still  unreachable 
at  the  terminal  even  after  the  battery  is  rested,  due  to  the  hys¬ 
teresis  voltage  h[k\. 

The  hysteresis,  also  known  as  “voltage  pull”  -  a  derivative 
expression,  depends  on  the  direction,  magnitude  and  duration  of 
the  preceding  current  as  well  as  on  the  SOC.  Hysteresis  is  known  to 
result  from  thermodynamical  entropic  effects,  mechanical  stress, 
and  microscopic  distortions  within  the  active  electrode  materials, 


which  perform  a  two-phase  transition  during  lithium  insertion / 
extraction  6].  Thus,  knowing  the  exact  value  of  hysteresis  is 
impractical.  However,  since  the  direction  of  the  hysteresis  voltage  is 
opposite  during  charging  and  discharging,  the  OCV  (at  a  certain 
SOC)  can  be  approximately  estimated  by  averaging  the  measured 
terminal  voltage  during  charging  and  discharging  when  the  entire 
experiment  is  performed  using  the  same  magnitude  of  current. 
Usually,  a  smaller  current1 2 3  is  required  in  order  to  reach  the  OCV  at 
both  ends,  i.e.,  to  measure  OCV  from  Vo(0)  to  VQ(1). 

The  SOC  of  the  battery  for  OCV-SOC  characterization  is 
computed  through  Coulomb  counting,  which  keeps  track  of  the 
amount  of  Coulombs  extracted  from/inserted  into  the  battery  and 
computes  the  SOC  as  a  ratio  of  remaining  Coulombs  and  battery 
capacity.  Hence,  the  knowledge  of  battery  capacity  is  crucial  for 
accurate  OCV-SOC  characterization.  Usually,  the  OCV-SOC  char¬ 
acterization  is  performed  on  a  new  battery,  hence,  manufacturer 
specified  capacity  (also  known  as  “rated  capacity”)  or  Coulomb 
counting  from  a  fully  charged  battery  to  a  fully  empty  battery  at 
very  low  current  values  will  serve  as  a  fairly  accurate  estimate  of 
the  initial  battery  capacity. 

However,  the  battery  capacity  is  known  to  fade  over  time; 
hence,  the  BFG  needs  to  account  for  the  capacity  fade  over  time.  At 
this  point,  the  OCV-SOC  characterization  made  using  incorrect 
capacity  will  become  more  and  more  unsuitable.  We  will  elaborate 
on  this  point  later  in  this  paper. 

The  conclusion  from  the  existing  literature  is  that  the  OCV-SOC 
characterization  changes  with  battery  temperature.  Hence,  many 
existing  approaches  suggested  computing  OCV  parameters  at 


1  It  must  be  noted  that  the  hysteresis  voltage  reversal  is  delayed,  i.e.,  when  the 
current  reversal  occurs,  the  hysteresis  reversal  follows  with  a  delay. 

2  Using  high  current  in  OCV  experiments  will  prevent  reaching  full  battery 
(SOC  =  1)  or  empty  battery  (SOC  =  0). 

3  A  BFG  needs  the  battery  capacity  information  in  the  “Coulomb  counting 
equation”  or  process  model,  see  Ref.  [3]. 
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Fig.  2.  Equivalent  circuit  of  a  battery  when  it  is  experiencing  dynamic  current  vs.  when  it  is  rested.  The  measurable  voltage  across  battery  terminals  is  indicated  by  v[/<]  and  the 
current  through  the  battery  is  /[/<].  The  currents  through  iq,  ...,  RN  are  indicated  by  i\[k],  ...,  ijv[/<],  h[k\  denotes  hysteresis  and  R0  denotes  the  series  resistance. 


intermediate  temperatures  spanning  possible  battery  usage.  The 
OCV  parameters  at  any  other  temperature  are  then  obtained 
through  interpolation.  In  these  characterizations,  the  SOC  is 
computed  relative  to  a  fixed  (nominal)  capacity  of  the  battery, 
which,  based  on  the  manufacturer  specifications,  holds  true  only  at 
the  “nominal  temperature”  of  the  battery. 

Hence,  one  of  the  focal  points  of  this  research,  backed  by  data 
collected  from  34  battery  cells  each  at  up  to  16  different  tempera¬ 
tures  spanning  -25  °C  to  50  °C  at  equal  steps,  is  to  investigate  if  the 
knowledge  of  battery  capacity  at  different  temperatures  plays  a  role 
in  OCV-SOC  characterization.  We  computed  the  OCV-SOC  char¬ 
acterization  curves  at  each  temperature  and  analyzed  the  differ¬ 
ences.  Further,  rather  than  using  the  rated  capacity  of  the  battery, 
we  computed  the  battery  capacity  through  slow  discharge  from  full 
to  empty  at  each  temperature.  The  experiments  were  repeated  on 
slightly  older  batteries  as  well.  Our  conclusion  is  that  the  OCV-SOC 
characterization  curve  remains  the  same,  regardless  of  the  tem¬ 
perature  or  age  of  the  battery,  as  long  as  the  correct  battery  ca¬ 
pacity  is  used  in  the  characterization.  The  implication  of  this 
conclusion  is  significant:  There  is  no  need  to  repeat  OCV-SOC 
characterization  experiments  at  different  temperatures— doing  so 
at  one  temperature  is  sufficient  for  SOC  tracking  at  all  temperatures 
of  battery  operation  as  long  as  the  capacity  is  estimated  dynami¬ 
cally  and  continuously  [1-3  . 

Once  the  pairs  of  {OCV,  SOC}  values  are  collected,  the  next  step 
is  to  model  these  data  given  later  in  (1 )  as  a  function.  From  a  mobile 
computing  perspective,  the  following  constraints  are  relevant  in 
selecting  such  a  function: 

•  Fewest  possible  parameters 

•  Ease  of  computation  and  inversion. 

Nowadays,  rechargeable  Li-ion  batteries  are  used  in  many  small 
electronic  devices  (e.g.,  smart  watches,  Google  glass,  etc.)  and  BFG 
is  a  desired  capability  in  these  devices,  however,  it  must  be  able  to 
function  with  limited  amount  of  memory  and  computational  needs. 
Since  many  of  the  BFG  algorithms,  such  as  the  ones  in  Refs.  [1—3], 
require  the  computation  of  OCV  and  its  derivative  for  a  certain  SOC 
and  finding  the  SOC  given  a  certain  OCV  (inversion  of  (1)),  the 
selected  OCV-SOC  characterization  function  needs  to  be  easily 
computable  and  invertible.  The  functions  that  we  eventually  select 
involve  logarithms  -  we  employ  a  computationally  efficient, 
numerically  stable  approach  to  implementing  logarithms  through  a 
Pade  approximation. 


4  Our  analysis  on  aged  batteries  was  done  using  a  limited  amount  of  data,  hence, 
the  conclusion  about  the  aged  batteries  must  be  considered  preliminary. 


The  novel  contributions  and  results  of  this  paper  can  be  sum¬ 
marized  as  follows: 

1.  Normalized  OCV  modeling.  We  show  that  the  OCV-SOC  charac¬ 
terization  is  the  same  for  any  temperature;  hence,  doing  so  at 
one  convenient  temperature  is  adequate  for  BFG  at  all  temper¬ 
atures.  We  also  show  some  preliminary  results  that,  when 
normalized  OCV  modeling  is  employed,  the  OCV-SOC  charac¬ 
terization  is  unaffected  by  aging. 

2.  New  OCV-SOC  models.  We  analyze  many  models  for  OCV-SOC 
characterization  and  propose  novel  models  that  are  suitable  for 
mobile  applications,  i.e.,  with  modest  memory  and  computa¬ 
tional  complexity. 

3.  Novel  approaches  for  OCV  model  parameter  estimation.  We  pro¬ 
pose  several  approaches  for  OCV  parameter  estimation. 

4.  Performance  comparison  of  OCV-SOC  modeling  as  it  relates  to  BFG. 
We  employ  several  performance  evaluation  metrics  in  order  to 
assess  the  suitability  of  models  in  BFG.  First,  the  required 
computational  complexity  of  each  model  is  summarized.  Then, 
OCV  modeling  accuracy  is  compared  at  a  certain  complexity. 

Some  recent  works  [8-10]  attempt  to  estimate  the  OCV-SOC 
parameters  online.  Online  estimation  of  the  OCV  parameters  helps 
to  develop  a  BFG  that  is  independent  of  characterization  and  hence 
applicable  to  any  type  of  battery.  A  challenge  in  online  OCV-SOC 
characterization  is  that  the  algorithms  in  Refs.  [8-10  rely  on  the 
end  usage  data  for  characterization;  continuous  data  between  full 
(SOC  =  1)  and  nearly  empty  (SOC  ~  0)  state  of  the  battery  is 
required  for  complete  characterization;  obtaining  these  data  will 
require  the  end  user  to  fully  charge  and  fully  drain  the  battery  at 
least  once;  partial  OCV-SOC  characterization  will  affect  the 
robustness  of  the  BFG  algorithm.  In  our  robust  BFG  algorithm 
(summarized  in  Fig.  1)  we  assume  that  the  OCV-SOC  parameters 
are  estimated  offline  and  the  equivalent  circuit  model  (ECM)  pa¬ 
rameters  and  the  capacity  are  estimated  online;  we  model  hyster¬ 
esis  as  a  bias  and  estimate  it  online  as  well.  The  battery  capacity, 
ECM  parameters  and  hysteresis  drift  with  environmental  and  usage 
patterns  and  age  of  the  battery;  estimating  these  parameters  re¬ 
quires  full  life  data  (which  is  typically  several  hundreds  of  cycles) 
under  different  environmental  and  usage  scenarios.  On  the  other 
hand,  the  OCV-SOC  parameter  estimation  requires  just  one  cycle  of 
data;  further,  we  show  in  this  paper  that  the  OCV-SOC  parameters 
are  stable  in  the  presence  of  temperature  changes  and  aging  of  the 
battery.  Hence,  offline  estimation  of  the  OCV  parameters  serves  our 
purpose  of  developing  a  robust  BFG  algorithm.  Further,  it  must  be 
noted  that  the  OCV  parameters  are  needed  for  the  (online)  esti¬ 
mation  of  ECM  parameters,  capacity  and  hysteresis;  accurate,  off¬ 
line  estimated  OCV  parameters  help  to  ensure  the  robust  online 
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estimation  of  these  quantities.  In  order  to  generalize  the  BFG  for 
multiple  battery  chemistries,  we  propose  in  Ref.  [11]  to  store  a 
carefully  selected  set  of  OCV  parameters  and  employ  probabilistic 
data  association  (PDA)  in  order  to  associate  the  correct  OCV 
parameter  to  the  battery  to  be  gauged. 

The  importance  of  the  total  capacity  of  the  battery  in  OCV-SOC 
parameter  estimation  is  not  fully  realized  in  the  literature.  The  data 
for  OCV-SOC  model  parameter  estimation  is  obtained  as  follows: 
the  OCV  data  is  collected  either  by  slow  discharge  or  by  intermit¬ 
tent  resting  12];  the  SOC  data  is  computed  through  Coulomb 
counting  which  requires  the  value  of  battery  capacity.  It  is  impor¬ 
tant  to  note  that  the  correct  value  of  battery  capacity  must  be  used 
in  OCV-SOC  modeling;  the  OCV-SOC  curve  will  be  different  for 
different  values  of  capacity  used  in  the  computation  of  SOC.  It  is 
observed  in  Refs.  [13,14]  that  when  normalized  by  “specific  cell 
capacity”,  the  OCV-SOC  curve  remains  the  same  for  the  same 
battery  chemistry;  this  observation  also  implies  that  the  OCV-SOC 
curve  remains  the  same  for  (albeit  slightly)  aged  batteries.  We  have 
confirmed  both  of  these  observations  using  34  battery  cells  of  11 
different  chemistries  at  eight  different  temperatures  and  two 
different  aging  conditions.  The  battery  capacity  changes  with 


Vo(s)=m  (i) 

where  VG(s)  denotes  the  OCV  and  s  <e  [0,1  ]  denotes  the  SOC.  In  this 
section,  we  summarize  important  OCV-SOC  functions. 

2.2.  Linear  regression-based  models 
A  linear  OCV  model  satisfies 

V0(s)  =  p(s)rk  (2) 

where  p (s)r  is  a  row  vector  of  linear/non-linear  functions  of  s  and  k 
is  the  OCV  parameter  vector.  It  must  be  noted  that  the  above 
expression  is  linear  in  the  parameter  vector  k,  not  necessarily  in  s. 

We  consider  the  following  linear-in-parameter  models  in  our 
analysis  with  the  objective  of  finding  an  appropriate  model  that 
satisfies  the  restrictions  in  the  number  of  parameters  and  compu¬ 
tational  complexity  while  providing  sufficient  modeling  accuracy 
[3,15-21]: 
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(i)  Shepherd  model  [15, 16] 

(ii)  Unnewehr  universal  model  [15,17] 

(iii)  Nernst  model  [15,18,  19] 

(iv)  Combined  model  [20] 

(v)  Combined  +  3  model  [3] 

(vi)  Polynomial  model  [21] 

(vii)  Exponential  model 
(viii)  Chebyshev  model 


(3) 


temperature;  when  normalized  by  the  capacity  at  the  temperature, 
the  OCV-SOC  curve  remains  the  same  across  temperature  changes 
as  well;  this  is  one  of  the  major  findings  in  this  paper.  We  also  detail 
how  the  normalization  should  be  done  by  total  capacity  of  the 
battery;  unlike  discharge  capacity,  the  total  capacity  is  the  same  at  a 
certain  temperature  at  any  discharge  rate. 

The  rest  of  the  paper  is  organized  as  follows:  A  summary  of 
existing  and  some  newly  proposed  OCV-SOC  models  are  described 
in  Section  2.  Section  3  discusses  the  different  OCV-SOC  model 
parameter  estimation  techniques.  The  novel  normalized  OCV 
modeling  is  summarized  in  Section  4.  In  Section  5,  we  elucidate  the 
different  metrics  for  the  selection  of  an  OCV-SOC  model  for  the 
BFG.  Section  6  evaluates  the  different  models  based  on  HIL  data 
collected  across  different  batteries  and  varying  operating  temper¬ 
atures.  Finally,  Section  7  concludes  with  a  summary  and  future 
research  directions. 

2.  OCV-SOC  models 

The  OCV  of  the  battery  has  a  monotonic  relationship  to  the  SOC. 
This  relationship,  usually  modeled  as  a  function,  is  the  backbone  of 
a  BFG.  In  a  general  form,  it  is  written  as 


k  =  [K0K,  K2...Km]t  (4) 

where,  M  +  1  is  the  length  of  the  parameter  vector  k  (for  example 
M  =  1  for  the  Shepherd  Model  and  M  =  L\  +  L2  for  the  polynomial 
model),  L  is  the  order  of  the  Chebyshev  polynomial  and  T\  is  the  2th 
Chebyshev  polynomial. 

Next  we  discuss  some  of  the  above  models. 

2.2.2.  Combined+3  model 

The  “Combined+3  model”,  first  used  by  the  authors  in  Refs. 
[  1  -3  ,  was  formed  by  including  three  additional  terms  1  /s2, 1  /s3  and 
1  Is4  to  the  existing  combined  model.  These  terms  are  introduced  to 
capture  the  sharp  decline  (or  “knee”)  in  the  OCV  as  the  SOC  reaches 
zero.  It  is  noticed  that,  such  sharp  knees  are  found  in  the  majority  of 
the  batteries  except  for  a  few  battery  chemistries.  For  batteries 
without  such  knee  shape  in  the  OCV-SOC  characterization,  the 
Combined  model  alone  will  be  adequate  to  accurately  capture  the 
OCV  curve. 

2.2.2.  Chebyshev  model 

The  Chebyshev  model  is  proposed  in  this  paper  due  to  its 
numerically  stable  characteristics  in  fixed  point  implementations 
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and  its  known  property  of  minimizing  the  maximum  error  (i.e., 
uniform  error  across  all  SOC  values)  22,23]. 

Let  us  define  sc  e  [-1,1]  as 


sc  =  2s  1 

(5) 

The  function  l/G(s)  can  be  approximated  in  terms  of  the  Che¬ 
byshev  polynomial  as 

Vo^X^TKSck) 

2=0 

(6) 

where 

Ko=}X>°(5ck) 

k= 1 

(7) 

Ki  = 

k=  1 

(8) 

and 

sck-cos  (  2l  \  k  —  1,2,  ...,L 

(9) 

T,(sck)  =  cos  i(2/<2L1)TC.  i<L 

(10) 

2.1.3.  Clenshaw  recurrence  for  the  evaluation  of  Chebyshev  model 
Once  the  coefficients  of  the  Chebyshev  expansion  are  known, 
then  V0(s)  for  any  given  s  <e  [0,1]  can  be  computed  using  the 
backward  recurrence  relations 

yL+ i=o  (li) 


Fig.  3.  Computing  In  through  Pade  approximation.  Top  plot  shows  ln(x)  computed  by 
Pade  approximation  (18),  Taylor  series  approximation  and  Matlab  (version  2011b).  The 
error  plots  in  the  bottom  are  relative  to  the  Matlab  version. 

This  numerically  stable  approach  to  the  evaluation  of  functions 
expressed  in  Chebyshev  series  is  known  as  Clenshaw's  recurrence. 

2.1.4.  Pade  approximation  of  natural  logarithm 

For  OCV-SOC  models  that  include  natural  logarithm,  there 
needs  to  be  an  approximate  way  to  implement  them  in  smaller 
micro  controllers  and  hardware  logic  circuits.  We  use  the  following 
expression,  which  gives  a  very  accurate  approximation  of  natural 
logarithm  through  Pade  approximation  [24] 


IfipadeO 


0.01 81 2x5  -  0.30555x4  -  1.30555x3  -  2x2  +x 
0.00396X5  -  0.11904X4  +  0.83333x3  -  2.22222x3  +  2.5x  -  1 


(17) 


yt  =  o 


(12)  The  lnpade  above  can  be  implemented  efficiently  using  Horner's 

method  as  follows: 


((((137x  -  2310)x  +  9870)x  -  15120)x  +  7560)x 
nPade(  -  x)  -  ((((30x  -  900)x  -  6300)x  -  16800)x  +  18900)x  -  7560 


(18) 


yL-i  =I<L-i  ~yL+ 1  +  2yi.sc 

(13) 

yi  =  i<\  -y3+2y2sc 

(14) 

yo=I<o  -y2  +2yiSc 

(15) 

Vo(s)  =  ( K0  -y2)  +yiSc  =y0  -yisc 

(16) 

In  Fig.  3,  the  lnpade  approximation  is  compared  with  In  imple¬ 
mented  using  5th  order  Taylor  series  approximation  and  the  In 
implementation  by  Matlab.  It  shows  that  the  Pade  approximation  is 
more  accurate  than  Taylor  series  approximation  of  the  same  order. 


5  Indeed,  higher  order  Taylor  series  approximations,  up  to  more  than  15th  order, 
were  tested  and  the  conclusion  is  that  the  Taylor  series  approximation  did  not 
improve  the  accuracy  with  respect  to  (18). 
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2.2.  Nonlinear  regression-based  models 

In  this  section,  we  review  three  non-linear  models  presented  in 
the  literature  and  present  a  new  non-linear  OCV-SOC  model  using 
a  rational  approximant. 

For  nonlinear  modeling,  we  designate 

V’o,n(S,*Co)=  V0(S)  (19) 

where 

K0  =  [kq  K\  . .  .km]  (20) 

is  the  nonlinear  OCV  parameter  vector  [25,27  . 


When  p*Q  (s)  and  pi  (s)  represent  two  straight  lines,  the  mini¬ 
mum  number  of  parameters  will  be  five:  two  parameters  each  for 
the  straight  lines  and 

2.4.  OCV-SOC  table 

An  alternate  approach  to  OCV  modeling  is  to  store  values  of  OCV 
and  SOC  as  a  table.  Table  1  shows  such  an  OCV-SOC  table  made  of 
11  SOC  points.  It  will  become  apparent  after  Section  4  that,  only  one 
OCV-SOC  table  is  all  that  is  required  for  fuel  gauging  at  all  tem¬ 
peratures.  Generally,  the  BFG  algorithm  requires  the  derivative  of 
the  OCV  function  to  be  used  in  an  extended  Kalman  filter-based 
SOC  tracking  approach;  hence,  the  table  needs  to  store  the  de¬ 
rivatives  as  well. 


Vo  ,n(S,  Ko) 


K0  +  /C-jS  -f-  K2  —  6  ^ 


_^5 

e  i-s 


K1  +K7e-K bC  ~s) 


Kq  ~  —  +  K2e 


/CqC  _|_  £ 2  -)-  /C3S  K/^S^  -\~  K$S^ 

E™  o  ^ 

1  +  E/=l  Kj+m ^ 


7,  m  >  0,n>0 


(a)  Double  exponential  model  [25] 

(b)  Nonlinear  exponential  model  1  [26] 

(c)  Nonlinear  exponential  model  2  [27] 

(e)  Rational  approximant 


(21) 


Rational  Approximant:  In  (21),  we  propose  a  model  based  on 
rational  approximant  which  gives  very  good  modeling  accuracy. 
The  model  based  on  rational  approximant  is  proposed  based  on 
Pade  approximation  ideas  of  24  .  The  coefficients  ko,  *1,  •••»  Km+n  of 
the  rational  approximant  can  be  uniquely  determined.  This  method 
is  also  known  as  m/n  order  rational  approximant. 


Table  1  was  formed  by  uniformly  sampling  SOC.  A  nonlinear 
sampling  can  be  employed  for  better  accuracy.  For  example, 
selecting  the  number  of  samples  in  a  “region”  proportional  to 
dV0(s)/ds  will  be  more  accurate  than  a  table  made  of  samples 
selected  linearly  across  s  e  [0,1  ]. 


2.3.  Hybrid  or  piecewise  linear  models 

It  was  noticed  that  using  a  single  function  to  model  the  entire 
shape  of  Vo(s)  requires  more  parameters  to  reach  a  certain  accuracy. 
The  same  level  of  accuracy  can  be  obtained  with  fewer  parameters 
through  hybrid  modeling.  The  objective  is  to  simplify  the  models  as 
much  as  possible,  hence,  only  linear  regression-based  models  are 
considered  in  this  approach. 

The  intuition  behind  hybrid  modeling  is  to  approximate  the  OCV 
models  as  a  few  linear  curves.  For  example,  for  some  batteries, 
starting  from  a  fully  charged  state,  the  OCV  decreases  gradually  (and 
linearly)  with  SOC  until  a  certain  (low)  SOC.  After  that,  the  OCV  drops 
steeply  (and  linearly)  with  decreasing  SOC.  This  particular  battery 
can  be  represented  using  two  straight  lines:  one  for  the  gradual  drop 
and  the  other  for  the  steep  drop.  Some  other  OCV-SOC  characteristic 
curves  could  be  approximated  by  three  (or  more)  different  models. 
We  consider  only  two  of  the  following  general  form: 


3.  OCV  parameter  estimation  techniques 

In  this  section,  we  describe  approaches  for  estimating  the  pa¬ 
rameters  of  the  models  described  in  Section  2.  For  accurate  estima¬ 
tion,  we  need  the  {VG(s),s}  pairs  spanning  s  <e  [0,1].  We  collect  this 
data  by  discharging  the  battery  from  full-to-empty  and  then  charging 
it  back  from  empty-to-full  with  a  very  low  current  (C/30  to  C/40). 

Fig.  4  shows  the  equivalent  circuit  of  a  battery  when  it  is 
charged/discharged  with  a  constant  rate.  First,  we  define  the  SOC  at 
a  given  time  as 


at  time  tk 


(23) 


Considering  the  voltage  measurement  errors,  the  measured 
voltage  is  written  as 


Zv[tk]  =  vfo]  +  n„[tfc] 


(24) 


V.(s)  =  /*sft  ifS-? 
1  p'.  (s)Tk{  if  s  <  ? 


(22) 


where  each  of  p\  (s)  and  p,G  (s)  denote  one  of  the  eight  different  OCV 
models  in  (3).  In  this  paper,  we  chose  the  combined  model  for  p'Q  (s) 
and  a  2nd  order  polynomial  model  for  p/Q  (s).  It  is  possible  to  have 
i  =  j  with  the  exception  that  the  corresponding  parameters  k\  and 
kJQ  respectively,  have  to  be  different. 


where  nv[tk]  is  the  voltage  measurement  noise  which  is  modeled  as 
white  Gaussian  with  standard  deviation  (s.d.)  av.  Since  the  current 
during  the  OCV  experiment  is  usually  kept  constant  through  highly 
accurate  programmable  charge/load  devices,  we  assume  that  there 
is  no  noise  involved  in  i[t/<]. 


6  C/N  is  a  common  representation  for  the  discharge/charge  rate  of  a  battery, 
where  C  denotes  the  battery  capacity  in  Ah  and  N  is  in  hours. 
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Table  1 

A  sample  OCV— SOC  table  that  requires  11  memory  space. 


s 

V.(s) 

dV0(s)/ds 

0.0 

3.0519 

44.4073 

0.1 

3.6594 

0.6830 

0.2 

3.7167 

0.8131 

0.3 

3.7611 

0.4807 

0.4 

3.7915 

0.4393 

0.5 

3.8275 

0.6055 

0.6 

3.8772 

0.8113 

0.7 

3.9401 

0.9777 

0.8 

4.0128 

1.0915 

0.9 

4.0923 

1.1822 

1.0 

4.1797 

1.3405 

During  the  OCV  experiment,  the  terminal  voltage  can  be  written 
as 

zv[tk]  =  V.  (sfoj)  +  h[tk]  +  i[tk\R0  +  nv[tk]  (25) 

where  is  the  hysteresis  or  voltage  “pull”  which  is  a  function  of 
current  and  SOC  of  the  battery.  Since  the  OCV  test  is  performed  at  a 
very  low  current,  we  assume  that  the  hysteresis  is  proportional  to 
the  current  only,  i.e. 


h[tk\*i[tk\ 

(26) 

Hence,  (25)  can  be  rewritten  as 

zv[tk]  =  V „  ( s[tk ])  +  i[tk\R0,h  +  nv[tk ] 

(27) 

where  the  effective  resistance 

>3 

O 

=r 

O 

+ 

(28) 

is  the  summation  of  the  battery  series  resistance  Ro  and  the  con¬ 
stant-current  hysteresis  equivalent  resistance ,  Rp. 


3.1.  Linear  model  parameter  estimation 

The  parameters  of  the  linear  OCV-SOC  model  in  (2)  can  be 
written  as, 


Fig.  4.  Equivalent  circuit  of  a  battery  when  it  experiences  constant  current.  When  the 
current  is  constant  for  a  long  time,  the  internal  capacitors  of  the  battery  becomes 
saturated  and  the  battery  equivalent  circuit  becomes  just  a  resistor. 


/ 

1 - 

4-1 

1 _ 

N 

p.  (s[tfc])T  i 

fc 

ko 

Ro,n  

+nv  t/< 

P  {ti<]T  k 

By  considering  a  batch  of  N  voltage  observations,  (29)  can  be 
written  as 

v  =  P/<  +  n  (30) 

where 

v=  [v[ti]  v[t2]...v[tN]]T  (31) 

p=  [p[ti]p[f2]...p[tN]]r  (32) 

n=  [n[fi]  n[t2]...n[tN]]r  (33) 

k  =  [k„  i?of!]ran  (M  +  2)vector  (34) 

The  least  squares  estimate  of  the  parameter  vector  is  given  by 

k=(pTp)~Vv  (35) 

Now,  for  a  given  SOC  s,  the  corresponding  OCV  estimate  V  Q  (s)  is 
computed  as 

V.{s)  =p„(s)7k  (36) 


where  k  is  formed  by  the  first  M  +  1  (4)  elements  of  k. 

In  Fig.  5,  the  estimated  values  of  /?o,h.  %  and  Rh  are  shown  using 
data  collected  on  a  Samsung  battery  with  serial  number  EB575152 
at  different  temperatures.  The  average  modeling  error  at  each 


Fig.  5.  Estimated  resistance  values.  Top:  the  estimated  value  of  effective  resistance 
Roh  =  R0  +  Rh  is  compared  with  the  series  resistance  R0,  which  was  computed  by 
applying  a  discharge  pulse  on  the  battery.  The  computed  value  of  hysteresis  equivalent 
resistance,  R j,  =  R0j 7  -  Ro,  is  also  plotted;  R j,  shows  a  similar  temperature  dependence 
as  R0.  Bottom:  corresponding  modeling  error  relative  to  the  voltage  observations. 
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temperature  is  also  shown  in  the  bottom  of  the  plot.  The  average 
modeling  error  is  discussed  later;  briefly,  the  average  modeling 
error  is  computed  by  averaging  the  voltage  error  for  the  entire 
range  of  SOC,  for  example,  the  error  at  50  °C  is  computed  by 
averaging  the  error  shown  in  the  bottom  Fig.  17.  The  box  plot  in 
Fig.  18  shows  this  modeling  error  across  multiple  batteries  and 
temperatures. 

3.2.  Total  least  squares  parameter  estimation 

In  contrast  to  the  classical  least  squares  (LS)  problem,  in 
which  the  solution  to  the  linear  system  of  equations  Pk  ~  v 
assumes  noise  only  in  the  observation  vector  v,  the  total  least 


squares  (TLS)  problem  assumes  that  the  elements  of  both  P  and  v 
are  contaminated  by  noise.  This  leads  to  the  following 
formulation 


kTLs  =  arg  min 

k,P,V 


(37) 


subject  to 


Pk  =  v  (38) 

where  [P  v]  and  [P  v]  are  matrices  of  equal  size.  The  matrix  P  is 
noisy  due  to  errors  in  computed  SOC,  as  discussed  in  Section  4. 
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Fig.  6.  Estimated  battery  capacities  at  different  temperatures.  The  capacities  are  computed  by  C/30  discharge.  Capacity  measurements  from  34  battery  cells  are  plotted  against 
temperature.  Markers  indicate  the  measurement  and  the  lines  connect  the  measurements  that  correspond  to  a  certain  cell.  Details  from  each  cells  are  shown  in  Figs.  7-10  for  clarity. 
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Fig.  7.  New  vs.  old  battery  capacity  against  temperature.  Only  cell  1  data  from  Fig.  6  is  shown. 
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Fig.  8.  New  vs.  old  battery  capacity  against  temperature.  Only  cell  2  data  from  Fig.  6  is  shown. 
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Fig.  9.  New  vs.  old  battery  capacity  against  temperature.  Only  cell  3  data  from  Fig.  6  is  shown. 
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Fig.  10.  New  vs.  old  battery  capacity  against  temperature.  Only  cell  4  data  from  Fig.  6  is  shown. 


Further  discussion  of  the  TLS  approach  and  its  method  of  solution 
can  be  found  in  Refs.  [28-30  . 


3.3.  Nonlinear  model  parameter  estimation 

For  nonlinear  models,  we  rewrite  (39)  in  the  following  form 
v(#c)  =  Vo  (ko  )  +  i  Roh  +  w  (39) 

where 

k  =  [/c.  J?0,h]T  (40) 


Fig.  11.  OCV-SOC  characterization  curves  of  new  vs.  aged  batteries.  New  battery 
curves  are  plotted  in  solid  blue  and  aged  battery  curves  are  plotted  in  dashed  red. 
Different  curves  of  the  same  type  correspond  to  temperatures  ranging  from  -25  °C  to 
50  °C.  The  new  and  aged  battery  OCV  curves  overlap.  (For  interpretation  of  the  ref¬ 
erences  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this 
article.) 


V.(Ko)  =  [\Zo,n(s[t1],K„)...VOi„(s[tN],K0]T  (41) 

i  =  [i[fi]  i'N •■■.%]]  (42) 

and  w  is  the  noise  vector. 

The  coefficients  of  the  nonlinear  regression-based  models  were 
computed  using  the  Matlab  optimization  toolbox  function  for 
nonlinear  least  squares  Isqnonlin.  The  nonlinear  least  squares 
problem  solves  the  following  problem 

k  =  argmin||v  -  v(#c)||  (43) 


Temperature  (°  C) 

Fig.  12.  A  quadratic  model  fitting  on  total  capacity  estimates:  Each  measurement  plot 
in  Fig.  6  is  fitted  with  a  quadratic  model.  Here  we  demonstrate  quadrature  model 
fitting  on  3  different  cells  of  the  same  chemistry  and  type.  The  corresponding  model 
capacity  is  plotted  in  dashed  line  of  the  same  color.  The  estimated  nominal  temper¬ 
ature  (60)  T nom  that  corresponds  to  the  maximum  capacity  is  shown  in  dash-dot  line 
of  the  same  color.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the 
reader  is  referred  to  the  web  version  of  this  article.) 
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Computed  nominal  temperature  (°  C) 


Fig.  13.  Analysis  of  estimated  nominal  temperature:  estimated  nominal  temperature 
of  34  battery  cells  of  different  manufacturers  is  used  for  this  experiment. 

3.4.  Hybrid  model  parameter  estimation 

It  is  well  known  that  OCV  is  a  monotonically  increasing  function 
of  SOC.  Consequently,  we  imposed  the  inequality  constraint  d VQ(s)/ 
ds  >  0  in  order  to  obtain  an  OCV  model  with  monotonic  estimates 
of  OCV  with  respect  to  SOC. 

Furthermore,  in  hybrid  models,  the  OCV  model  is  assumed  to 
follow  a  particular  function  until  s  =  £  and  then  assumes  a  different 
model  structure.  In  order  to  realize  such  a  model  structure,  the  least 
squares  problem  is  implemented  with  additional  equality  con¬ 
straints  as  shown  below: 


ki  -  arg  min 

ki 


S<E[0,?] 


(44) 


subject  to 


Fig.  14.  Normalized  OCV  modeling.  Measured  terminal  voltage  is  shown  in  solid  and 
computed  OCV  is  shown  in  dashed  lines.  The  solid  lines  show  the  measured  voltage 
and  the  dashed  lines  show  the  normalized  OCV  curves;  the  normalized  OCV  curves 
nearly  coincide  for  all  temperatures. 


Fig.  15.  Nominal  vs.  normalized  OCV  modeling.  Top:  nominal  OCV  modeling  uses 
Cn0m  =  1.5  Ah  in  computing  SOC  at  all  temperatures.  Bottom:  normalized  OCV 
modeling  computes  Cbatt  at  each  temperature  before  computing  SOC. 


dp*k 

ds 


>0 


(45) 


Pi(s)k1  _  -  p'(s)k2 


s=C 


s=C 


0 


(46) 


dp'fsjk 


ds 


dp'(s)  k: 


ds 


0 


(47) 


dV(s)k- 


ds2 


dV(s)k2 


ds2 


0 


(48) 
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dV(s)k2 


s=C 


ds2 


0 


(53) 


The  constrained  least  squares  algorithm  software  in  the  opti¬ 
mization  toolbox  of  Matlab  is  used  for  the  optimizations  above. 


4.  Normalized  OCV  modeling 


So  far  in  OCV  model  parameter  estimation  the  SOC  is  assumed  to 
be  known,  for  example,  the  SOC  in  (27)  is  computed  through 
Coulomb  counting: 


k 


k 


V^k  fifa-i] +ifa] 

3600C V  2 


(54) 


where 


Fig.  16.  OCV-SOC  characterization  data  and  sample  OCV  plot.  Top  two  plots  indicate 
the  applied  current  i[t/J  and  measured  voltage  zv[t/<].  The  plot  at  the  bottom  shows  the 
same  zv[t/< ]  in  the  middle  plot  above  against  SOC,  instead  of  time.  The  SOC  was 
computed  by  Coulomb  counting,  explained  later  in  (54).  The  OCV  is  computed  using 

(36). 


k2  -  arg  min 
1<2 


p'(s)k2-v)  se(?,  1 


(49) 


subject  to 
dp,k2 


ds 


>0 


(50) 


P*(s)ki  -p'(s)k2 


s=C 


0 


(51) 


dp'fsjk 


ds 


dp7  (s)k: 


ds 


0 


s=C 


(52) 


—  k  k— i 

is  the  sampling  interval,  and 

_  {  v c  <  1  during  charging 
71  ~  |  nd  =  1  during  discharging 

is  the  efficiency  and  C  denotes  the  battery  capacity  (in  Ampere- 
hours),  which  plays  a  major  role  in  OCV-SOC  characterization  of 
the  battery,  as  we  describe  below. 

There  are  two  different  ways  of  characterizing  OCV  based  on  the 
assumption  of  the  battery  capacity: 

r  _  f  ^nom  Nominal  OCV  modeling  fl-7) 

\  C5att  Normalized  OCV  modeling  ^ 

where  Cbatt  is  the  total  battery  capacity  (clearly  defined  in  a  sub¬ 
sequent  paragraph). 

We  define 

1.  Nominal  OCV  modeling ,  when  the  same  battery  capacity  C  =  Cn0m 
is  used  at  all  temperatures  for  computing  the  SOC  in  (27)  and 

2.  Normalized  OCV  modeling ,  when  the  computed  total  battery 
capacity  C  =  Cbatt  at  each  temperature  is  used  for  computing  the 
SOC  in  (27). 

In  Fig.  6,  we  show  the  total  capacities,  computed  by  C/30 
discharge  from  full  to  empty,  of  cells  listed  in  Table  2  at  different 
temperatures.  All  the  batteries  are  purchased  through  commercial 
vendors;  the  manufacturer  data  sheets  were  not  available  to  any  of 
the  batteries.  The  chemical  composition  of  the  batteries  was 
learned  through  destructive  analysis  on  sample  batteries.  All  the 
Samsung  cells  were  found  to  be  of  LiCo02  cathode,  Carbon  anode 
systems,  with  carbonate-based  electrolytes  containing  LiPF6  as  the 
primary  salt.  The  capacity  vs.  temperature  is  observed  to  take  a 
concave  form  in  general,  i.e.,  the  maximum  total  capacity  is 
achievable  only  at  a  particular  temperature. 

Now,  let  us  define  the  following  terms  related  to  battery 
capacity: 

•  Maximum/full-charge  battery  open  circuit  voltage  (OCVmax): 
OCVmax  refers  to  the  OCV  of  the  battery  when  it  is  fully  charged, 


(55) 


(56) 


B.  Pattipati  et  al.  /  Journal  of  Power  Sources  269  (2014)  317—333 


329 


Temperatuer  =  50°C 


Fig.  17.  Sample  OCV  modeling  error.  Top:  Measured  voltage  and  modeled  OCV.  Bottom: 
OCV  modeling  error. 


OCVmax  =  Vo  (s)|s=1  (58) 

It  must  be  noted  that  due  to  the  internal  impedance  of  the 
battery  and  hysteresis,  the  measured  cell  voltage  of  the  battery  will 
be  higher  than  OCVmax  right  after  the  battery  is  fully  charged. 

•  Minimum/empty  battery  open  circuit  voltage  (OCVmin): 

Similar  to  (58) 

OCVmin  =  VD  (s)|s=0  (59) 

where,  due  to  the  internal  impedance  of  the  battery  and  hysteresis, 
the  measured  cell  voltage  of  the  battery  will  be  lower  than  OCVmin 
right  after  the  battery  is  fully  discharged. 

•  Rated  battery  capacity,  Cr:  Rated  capacity  is  the  manufacturer 
specified  (labeled)  capacity  of  the  battery.  Manufacturer  speci¬ 
fied  capacity  should  be  treated  as  an  approximate  measure  of 
the  true  value  of  battery  capacity  at  the  specified  temperature. 

•  Total  capacity,  Cbatt  or  Cbatt (T):  Total  battery  capacity  is  the 
maximum  amount  of  Coulombs  that  can  be  discharged  from  a 
battery  at  cell  temperature  T,  starting  from  a  fully  charged  bat¬ 
tery  corresponding  to  its  OCV  =  OCVmax  until  the  battery  is  fully 
discharged,  i.e.,  until  the  battery  terminal  reaches  OCV  =  OCVmin 
using  an  infinitesimal  load  current. 

•  Nominal  capacity,  Cn0rm  Nominal  capacity  is  the  maximum 
amount  of  Coulombs  that  can  be  discharged  from  a  battery  at 
any  cell  temperature.  In  other  words,  nominal  capacity  is  the 


7  The  series  resistance  of  the  battery  prevents  measured  terminal  voltage  from 
reaching  OCVmjn  and  hence  fully  discharging  the  battery.  A  smaller  load  current 
makes  sure  the  battery  is  fully  emptied  as  much  as  possible.  The  gradient  of  the 
OCV— SOC  curve  is  very  high  (see  Table  1  and  a  sample  OCV— SOC  plot  in  Fig.  11) 
close  to  an  empty  battery,  which  helps  minimize  the  errors  in  total  capacity  esti¬ 
mation  by  slow  discharging. 

8  Most  of  the  battery  literature  treats  nominal  capacity  as  exactly  similar  to  the 
rated  capacity.  We  introduce  this  slightly  altered  definition  based  on  observations 
in  the  literature  where  the  nominal  capacity  is  treated  as  the  maximum  possible 
capacity.  For  example,  in  Ref.  [31  ,  the  nominal  capacity  of  a  power  plant  is  the 
maximum  possible  power  output.  Nominal  capacity  is  technically  achievable,  but, 
not  always  practical.  In  the  battery  example,  the  nominal  capacity  is  only  achievable 
at  the  nominal  temperature. 


maximum  total  capacity  of  the  battery.  The  temperature  at 

which  the  total  battery  capacity  is  maximum  is  defined  as  the 

nominal  temperature  of  the  battery,  i.e., 

T  nom  =  arg  max  Cbatt(7)  (60) 

Cnom  —  Qoau  ( 7  nom)  (61) 

The  total  capacity  Cbatt  increases  with  temperature  up  to  the 
nominal  temperature.  Beyond  that,  the  total  capacity  starts  to 
decline  with  increasing  temperature.  The  above  definition  implies 
that  the  nominal  temperature  of  the  battery  varies  from  battery  to 
battery.  It  must  be  noted  that  the  manufacturer  specified  nominal 
temperature  (and  nominal  capacity)  are  only  approximate.  Hence, 
for  accurate  BFG,  the  nominal  temperature  of  the  battery  must  be 
dynamically  estimated  —  we  will  present  those  details  in  a  future 
publication. 

Table  3  lists  the  differences  between  nominal  and  normalized 
OCV  modeling. 

The  conclusion  is  that  the  normalized  OCV  modeling  has  many 
advantages  over  nominal  OCV  modeling. 

Fig.  12  shows  the  quadratic  model  fitting  of  the  total  capacity 
estimates  against  temperature  using  data  collected  from  3  different 
cells  of  the  same  chemistry  and  type  of  battery.  Once  the  model  is 
fitted,  the  nominal  temperature  and  nominal  capacity  of  the  battery 
can  be  estimated  based  on  (60)  and  (61). 

The  nominal  temperatures  are  estimated  for  all  the  cells  listed 
in  Table  2  and  Fig.  13  summaries  the  estimated  nominal  tem¬ 
peratures  as  a  histogram.  Based  on  these  (limited  number  of) 
observations,  the  most  cells  had  a  nominal  temperature  of 
22.5  °C  and  the  average  nominal  temperature  is  23.5  °C.  It  must 
be  noted  that  in  most  of  the  batteries  in  Table  2,  the  manufac¬ 
turer  specified  nominal  temperature  is  25  °C.  Hence,  for  superior 
accuracy,  a  BFG  must  estimate  (and  track)  the  Capacity  of  the 
battery. 

Fig.  14  shows  an  example  of  normalized  OCV  modeling  at 
different  temperatures.  It  can  be  seen  that  regardless  of  the 
measured  terminal  voltages  at  different  temperatures,  the  corre¬ 
sponding  OCV  curves  are  nearly  identical.  This  is  further  illustrated 
later,  using  Fig.  15. 

Fig.  15  shows  the  comparison  of  OCV-SOC  curves  obtained 
from  nominal  (top)  and  normalized  (bottom)  OCV-SOC  charac¬ 
terization.  The  nominal  OCV  modeling  used  Cnom  =  1.5  Ah  in 
computing  SOC  at  all  temperatures  and  the  normalized  modeling 
used  the  computed  total  capacity  at  each  temperature  in  order  to 
compute  SOC.  It  can  be  observed  that  in  normalized  OCV 
modeling,  the  curves  computed  at  all  temperatures  are  nearly 
identical. 

Fig.  11  shows  the  OCV-SOC  characteristic  curves,  obtained  on  a 
new  as  well  as  aged  battery  at  different  temperatures.  First,  the 
OCV  data  was  collected  from  a  newly  bought  battery 
at  -25  °C,  -15  °C,  -5  °C,  15  °C,  25  °C  and  35  °C  and  45  °C.  Then, 
the  battery  was  left  for  “self  aging”  for  approximately  one  year. 
Then,  the  OCV  data  was  collected  from  this  aged  battery 
at  -20  °C,  -10  °C,  0  °C,  10  °C,  20  °C  and  30  °C,  40  °C  and  50  °C.  It 
can  be  seen  that  no  significant  difference  of  patterns  was  found 
with  respect  to  the  age  and  temperature  of  the  battery.  Conse¬ 
quently,  a  single  normalized  OCV-SOC  curve  in  conjunction  with 
an  adaptive  battery  capacity  estimation  algorithm  can  be  used  in 


9  There  was  no  particular  reasons  for  the  selection  of  these  temperatures.  At  first, 
it  was  decided  to  do  the  tests  in  10  °C  intervals  and  later  it  was  decided  to  the  tests 
(on  aged  batteries)  at  slightly  different  temperature. 
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(a)  Combined+3  (b)  Polynomial  (c)  Exponential 


Double  Exponential  Model  Nonlinear  Exponential  Model  2 


Temperature  Temperature  Temperature 


(d)  Chebyshev 


(e)  Double  Exponential 


(f)  Nonlinear  Exponential  2 
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Fig.  18.  RMS  error  across  different  battery  cells.  RMS  value  of  the  OCV  modeling  error  across  different  battery  cells  is  given  as  a  box  plot  (the  central  (red)  mark  of  the  box  is  median; 
the  edges  of  the  box  are  25th  and  75th  percentile  values;  the  whiskers  extend  to  the  most  extreme  data  points  not  considered  outliers,  and  outliers  are  plotted  individually  with  a 
*+’)  at  each  temperature.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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Table  2 

Hardware-in-the-loop  battery  OCV— SOC  characterization  data  overview. 


Battery  Supplier 

Battery  type 

Number  of  batteries 

Rated  capacity  (mAh) 

Samsung 

EB575152VA 

4 

1500 

Samsung 

EBL1A2GBA 

4 

1650 

Samsung 

EB555157VA 

2 

1750 

Samsung 

EB504465VA 

4 

1500 

Samsung 

AB463651 

4 

960 

Samsung 

LI  G6LLAGS 

2 

2100 

Nokia 

BP-4L 

4 

1500 

LG 

LGIP-530B 

2 

1100 

RIM 

F-Sl 

2 

1270 

RIM 

M-Sl 

4 

1500 

Huawei 

HB4Q1 

2 

1670 

the  BFG  to  obtain  accurate  SOC  estimates  for  varying  usage  con¬ 
ditions  over  the  lifetime  of  a  battery. 


5.2.4.  Root-mean  square  (RMS)  error 


RMS  - 


l|v-v|| 

VN-M 


or 


Vmse 


(65) 


where  N  is  the  number  of  data  points,  M  is  the  number  of  param¬ 
eters,  v  is  the  predicted  value  of  v  using  the  estimated  parameters, 
for  example,  for  linear  models 


v  =  Pk 


(66) 


and 


1  N 

v=nE?(  '■) 

i=i 


(67) 


5.  OCV-SOC  model  metrics  for  selection  of  the  best  model  5.2.  Model  evaluation  metrics 


In  this  section,  we  develop  a  systematic  approach  to  select  the 
best  approach  for  OCV  models.  This  is  done  by  computing  several 
error  metrics  and  model  evaluation  metrics ,  described  below,  and 
then  by  ranking  the  models  based  on  Borda  counts. 

5.2.  OCV  prediction  error  metrics 

The  following  four  error  metrics  are  used  for  predicting  v, 
defined  in  (31)  to  evaluate  the  models: 

5.2.2.  Best  fit 


Model  evaluation  metrics  consider  the  trade  off  between  the 
number  of  model  parameters  and  the  number  of  data  points.  The 
following  four  metrics  are  important  ones. 

5.2.2.  Akaike's  information  criterion  (AIC) 

If  the  models  are  fitted  using  Least  Squares,  then  [32]  suggests 
the  following  analog  of  AIC: 

AIC  =  Nln  (^)  +  2(M  +  1 )  (68) 


BF(%) 


v  —  v||\ 

- =r  x  100 

v  ~  vll / 


5.2.2.  R2  fit 


x  100 


x  100 


-  [BF(2  —  BF)]  x  100 


5.2.3.  Max  error 


ME  =  max{|Vj  -  Vj|} 

l 


where 


(62) 


N 


i=l 


(69) 


with 

e  =  v  -  v  (70) 


In  the  above  S2  is  the  sum  of  the  squares  of  errors  (SSE),  e*  is  the 
tth  element  of  the  residual  vector  e  and  M  is  the  number  of  pa¬ 
rameters  in  the  OCV  model.  The  better  the  model,  the  lower  the  AIC. 


5.2.2.  Akaike’s  information  criterion  2  (AIC2) 

A  second  version  of  AIC,  given  below,  is  useful  when  N  »  M 


(64) 


AIC2 


In 


(71) 


Table  3 

Comparison  of  nominal  vs.  normalized  OCV  modeling. 


Nominal  OCV  modeling  Normalized  OCV  modeling 


The  SOC  at  any  temperature  is 
relative  to  the  fixed,  nominal 
capacity  of  the  battery 
There  is  no  need  to  know  the  total 
capacity  at  a  temperature  in  order 
to  compute  the  SOC 
Gives  accurate  (durable)  OCV 

parameters  only  with  a  new  battery 

Capacity  estimation  becomes 
erroneous  with  an  aged  battery 


The  SOC  is  relative  to  the  total 
capacity  of  the  battery  at  the 
given  temperature 
Need  to  know  the  total  capacity 
in  order  to  compute  SOC  at  a 
particular  temperature 
OCV— SOC  characterization  can 
be  done  using  an  aged  battery 
as  well 

Capacity  estimation  works  well 
even  for  an  aged  battery 


where  the  loss  function  is  defined  as 


eTe 

~N~ 


5.2.3.  Akaike’s  final  prediction  error  (FPE) 


FPE 


1  +  M 

1  +  N 

1  _  M 

1  N 


(72) 


(73) 
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5.2.4.  Bayesian  information  criterion  (BIC) 

The  derivation  of  BIC  assumes  equal  priors  on  each  model  and 
noninformative  priors  on  the  parameters,  given  each  model.  The 
goal  of  the  BIC  is  to  find  the  best  (i.e.  highest  posterior  probability) 
model  for  prediction. 


BIC  =  2(Ljv)  +  (M  +  l)lnN 


(74) 


The  negative  log-likelihood  given  the  pdf  of  the  residuals 
(assuming  normal  or  Gaussian)  conditioned  on  the  parameters  k 
and  the  s.d.  of  residuals  a  is  given  by 


Ln  =  -ln{L(k;  e)} 


+  0.51n(  2tt(j 


(75) 


where,  Ln  is  the  negative  log-likelihood,  L  is  the  likelihood,  k  is  the 
parameter  vector  which  minimizes  Ln ,  and  a  is  the  s.d.  of  the  re¬ 
siduals  e. 


5.2.5.  Minimum  description  length  (MDL) 


MDL  =  SBf 


,  Min  AT 

1 


(76) 


5.3.  Combining  rankings 

The  Borda  count  [33]  is  an  easy,  intuitively  appealing,  and 
powerful  method  for  combining  different  rankings.  In  this  paper, 
we  use  Borda's  count  method  to  combine  the  rankings  of  different 
model  selection  metrics.  The  Borda  count  is  originally  a  voting 
method  in  which  each  voter  gives  a  complete  ranking  of  all  possible 
alternatives.  The  highest  ranked  alternative  or  OCV  Model  (in  an  n- 
way  vote)  gets  n  votes  (or  points)  and  each  subsequent  alternative 
gets  one  less  vote  (i.e.  number  2  gets  n  -  1  votes  and  number  3  gets 
n  -  2  votes  and  so  on).  Then,  for  each  alternative,  all  the  votes  (or 
points)  are  added  up  and  the  alternative  (i.e.  OCV  Model)  with  the 
highest  number  of  votes  wins  the  election.  The  OCV  Models  are 
then  ranked  in  descending  order  of  votes  (or  points)  to  give  the 
overall  Borda  ranking. 

5.4.  Computational  complexity  metrics 

Most  of  the  present  day  applications  of  BFG  need  to  operate  with 
limited  amount  of  memory  and  computations.  This  metric  aims  to 
quantify  the 

1.  Memory:  number  of  variables  required  to  define  the  OCV-SOC 
characterization 

2.  Computations:  the  computational  complexity  in  finding  the 
OCV  for  a  specific  SOC. 

In  Table  4,  we  list  the  number  of  parameters  to  be  stored  in 
memory  as  well  as  the  computational  requirements  for  calculating 
OCV  for  a  given  SOC  in  terms  of  the  number  of  additions  (add), 
number  of  floating  point  multiplications  (mult)  and  the  number  of 
logarithmic  computations  (ln). 

6.  HIL  data  analysis 

The  OCV-SOC  characterization  data  was  gathered  by  hardware- 
in-the-loop  experiments  conducted  on  9  different  batteries  at  16 
different  temperatures  from  -25  °C  to  +50  °C  in  steps  of  5  °C  using 
a  Tenney  thermal  chamber  and  an  Arbin  BT-2000  cell  test  station. 
The  battery  manufacturers,  battery  types,  number  of  batteries  used 


Table  4 

Memory  and  computational  requirement  of  OCV-SOC  models. 


Model 

Required  memory 

Computation:  SOC  — >  OCV 

Shepherd  model  (2.1) 

2 

1  add  +  1  mult 

Unnewehr  model  (2.1) 

2 

1  add  +  1  mult 

Nernst  model  (2.1) 

3 

2  add  +  2  mult  +  21n 

Combined  model  (2.1) 

5 

4  add  +  4  mult  +  21n 

Combined+3  model  (2.1) 

8 

7  add  +  7  mult  +  21n 

Polynomial  model  (2.1) 

1+  1 

L  add  +  L  mult 

Exponential  model  (2.1) 

1+  1 

L  add  +  L  mult  +  Lexp 

Chebyshev  model  (2.1) 

L 

1  add  +  2  mult 

Double  exponential  (2.2) 

5 

5  add  +  3  mult  +  2exp 

Nonlinear  exponential  2  (2.2) 

6 

1  add  +  2  mult 

Rational  approximant  (2.2) 

2 

1  add  +  2  mult 

Hybrid  models  (2.3) 

Variable 

Model  dependent 

in  the  experiment,  and  their  rated  capacities  are  summarized  in 

Table  2. 

For  each  cell,  the  data  is  collected  in  the  following  manner:  First, 
the  cell  is  rested  at  the  test  temperature  at  full  charge.  Then,  the  cell 
was  discharged  at  C/30  rate  until  fully  discharged  to  the  minimum 
voltage  (3  V)  followed  by  a  short  rest  period.  Subsequently,  the  cell 
was  charged  at  C/30  rate  until  fully  charged  (4.2  V).  The  low  rates 
were  used  to  minimize  the  cell  dynamics.  Some  optional  charge/ 
discharge  pulses  (approximately  C/3  rate)  were  used  to  determine 
(initialize)  the  cell  resistances  and  time-constants.  The  OCV  was 
computed  by  averaging  the  cell  voltage  as  a  function  of  SOC  under 
discharge  and  charge.  This  eliminates  the  presence  of  hysteresis 
and  ohmic  resistance  effects  in  the  final  OCV  function. 

Fig.  16  shows  the  typical  OCV-SOC  characterization  data  that  is 
collected  from  one  of  the  cells  and  the  resulting  OCV  curve 
obtained. 

Fig.  17  shows  a  sample  OCV  model  obtained  for  a  particular 
battery  cell  (cell  4  of  Samsung  EB575152)  at  50  °C.  The  modeling 
error  at  a  particular  SOC  is  computed  by  subtracting  modeled  OCV 
from  the  average  of  measured  voltages  (during  charging  and 
discharging). 

In  Fig.  18,  we  give  a  comprehensive  summary  of  OCV  modeling 
showing  the  errors  of  different  OCV  modeling  approaches,  across 
different  batteries  and  temperatures.  First,  the  average  value  of  the 
OCV  modeling  error  of  an  each  battery  cell  at  each  temperature  is 
computed  by  averaging  the  OCV  modeling  error  across  all  SOC 
values  (i.e.,  by  averaging  the  OCV  modeling  error  shown  in  the 
bottom  plot  of  Fig.  17).  Then,  in  Fig.  18,  the  average  error  across  all 
the  battery  cells  listed  in  Table  2  is  shown  as  a  box  plot  for  each 
temperature.  The  Chebyshev  model  must  be  noticed  for  its  smallest 
error  across  all  temperatures  and  for  its  consistency  across  different 
batteries  as  indicated  by  the  (compact)  sizes  of  boxes  and  the  lack  of 
outliers  across  all  temperatures. 

The  results  of  the  model  evaluation  metrics  and  error  metrics 
have  been  combined  and  presented  in  fable  5.  The  Chebyshev 
model  performs  the  best  among  all  the  models  based  on  the  met¬ 
rics  presented  in  the  paper. 

7.  Conclusions 

This  paper  presented  several  approaches  and  analysis  related  to 
the  open  circuit  voltage  characterization  of  Li-ion  batteries  with  a 
special  emphasis  on  mobile  and  wearable  electronic  device 


10  Using  box  plot  function  in  Statistics  Toolbox  of  Matlab.  On  each  box,  the  central 
mark  is  the  median,  the  edges  of  the  box  are  the  25th  and  75th  percentiles,  the 
whiskers  extend  to  the  most  extreme  data  points  not  considered  outliers,  and 
outliers  are  plotted  individually. 
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Table  5 

Model  selection  metrics  rankings. 
OCV  Model  Rankings 


AIC  AIC2  FPE  BIC  MDL  RMSE  R 2  BF  Max  error  Borda 


Chebyshev  model 

1 

1 

1 

1 

1 

1 

1 

1 

4 

1 

Rational  approximant 

2 

2 

2 

2 

2 

2 

2 

2 

5 

2 

Polynomial  (4th  order) 

3 

3 

4 

3 

4 

3 

4 

3 

1 

3 

Hybrid  model 

5 

5 

3 

4 

3 

4 

3 

4 

2 

4 

Combined  model 

6 

6 

5 

6 

5 

5 

5 

5 

7 

5 

Exponential  model 

7 

7 

6 

5 

6 

6 

6 

6 

6 

6 

Nonlinear  exponential 

8 

8 

7 

7 

7 

7 

7 

8 

9 

7 

model  2 

Combined+3  model 

4 

4 

12 

11 

12 

7 

11 

7 

3 

8 

Nernst  model 

9 

10 

8 

8 

8 

9 

8 

9 

10 

9 

Unnewehr  universal 

11 

11 

9 

10 

9 

10 

9 

10 

11 

10 

model 

Double  exponential 

10 

9 

10 

9 

10 

11 

12 

12 

12 

11 

model 

Shepherd  model 

12 

12 

11 

12 

11 

12 

10 

11 

8 

12 

applications.  Memory  and  computational  power  becomes  scarce 
(due  to  prioritization)  for  BFG  in  these  applications,  yet,  the  BFG 
capabilities  such  as  prediction  of  time  to  shut  down  and  remaining 
useful  life  are  considered  desirable  features  to  have. 

Offline  characterization  of  open  circuit  voltage  is  crucial  for  BFG, 
especially  in  memory  and  computing  power  constrained  electronic 
applications.  We  include  several  existing  OCV  modeling  functions 
and  many  newly  proposed  OCV  modeling  approaches  in  a  perfor¬ 
mance  analysis  in  terms  of  OCV  modeling  accuracy,  computational 
complexity  and  memory  requirements. 

It  was  discovered  that  significant  savings  in  memory  can  be 
achieved  by  the  proposed  normalized  OCV  modeling,  in  which  the 
SOC  at  each  temperature  is  defined  in  terms  of  the  temperature 
specific  capacity.  As  a  result,  we  found  that,  the  OCV-SOC  charac¬ 
terization  remains  practically  the  same  regardless  of  the  tempera¬ 
ture  or  age  of  the  battery.  This  finding  is  crucial  for  online  battery 
capacity  estimation,  i.e.,  the  normalized  OCV  parameters  will  serve 
as  a  means  for  accurate,  online  battery  capacity  estimation,  such  as 
the  one  presented  in  Ref.  [2  . 

Apart  from  the  benefits  in  mobile  applications,  the  normalized 
OCV  modeling  approach  reported  in  this  paper  has  impact  on  other 
applications  such  as  in  electric  vehicles  and  energy  storage  grid 
where  accurate  capacity  estimation  is  critical. 
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